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1. Introduction 

Flavor physics is one of the urgent applications of Lattice QCD. However, the fact that the 
heavy quark masses are large in lattice units is a long-standing problem for heavy quark physics 
with LQCD. In the application, ma <C 1 is no longer true and the terms containing (ma) n (with a 
the lattice spacing) become significant. As a direct simulation with a <C 1/m costs too much, we 
resort to effective field theories. Various heavy quark effective actions were developed and used for 
different physical systems, see Refs. [^j, S BJ BJ| for reviews on this topic. 

In this proceeding, our work is based on the so called relativistic heavy quark (RHQ) action [||, 
|6[ [7|, ||]. The lattice form of the action, following the formulation proposed in [8j|, can be written 
as: 

5 = E^»| m o + 7oA)-^a^ + C y-b- X -a{bf - a£ ^c P a MV /> v j % (1.1) 

In the heavy quark case, the temporal covariant derivative Do is around the order of ma and should 
not be treated the same way as the spatial derivatives D,, which are of order Agcz)« or a s ma de- 
pending on the system under investigation. Following the Symanzik improvement procedure, we 
found that only the three free parameters mo, cp and £ need to be tuned to remove all errors of order 
(ma) n and \pa\. Thus, if the parameters are correctly tuned, the action will have small cutoff ef- 
fects: {Kqcdo) 2 for heavy-light systems and (a s ma) 2 for heavy quarkonium. The main purpose of 
this work is to determine the three parameters by matching to physical quantities for charmed sys- 
tems, making more accurate predictions for charmed mesons possible. The lattice spacing can also 
be obtained with reasonable precision if we treat it as a fourth quantity to be adjusted to correctly 
predict the charmed meson spectrum. All work has been done on dynamical 2+1 flavor lattices, 
which is a continuation of work done by H.-W. Lin |§J. 

The lattices used in this work are the dynamical 2+1 flavor 24 3 x 64 DWF lattice configurations 



generated by the RBC and UKQCD collaborations QlOp . For each configuration, we place sources 
at times 0, 16, 32 and 48 separately for better statistics; see Fig. [[} Part of the data was collected 
and the analysis was done during and after the lattice conference. This additional data is included 
in this proceeding for completeness. Binning the data every two configurations had no effects on 
the results, which suggests the auto-correlation of the lattice configurations is negligible. 



volume 


L s 


(m sea ,m s ) 


Traj(step) 


# of configs 


24 3 x 64 


16 


(0.005,0.04) 


900-4500(40) 


91 


24 3 x 64 


16 


(0.01,0.04) 


900-4500(40) 


91 


24 3 x 64 


16 


(0.02,0.04) 


1885-3605(20) 


87 


24 3 x 64 


16 


(0.03,0.04) 


1000-3060(20) 


104 



2. Determine the RHQ action and the lattice spacing 

To determine the action in such a way that errors are controllable, we tune the parameters 
by matching physical observables sensitive to them to their experimental values. The parameters 
are then determined for each ensemble with different light sea quark masses and extrapolated to 
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the chiral limit. The physical on-shell quantities we are going to use are mass combinations of 
pseudo-scalar (PS), vector (V), scalar (S) and axial-vector (AV) mesons in heavy-heavy (hh) and 
heavy-light (hi) systems [JTTJ] . 

• spin-averaged: mf a = \{mf s + ?>mf), mf a = \{m% + 2>m^) 

• hypernne splitting: m 1 ^ = m v h — rr&, m¥ = niy — m h p l s 

• mass ratio: where E 2 = m\ + j^_P 2 > m i'- rest mass, m^. kinetic mass. 

• spin-orbit averaged and splitting: m^ s = mf v — nif, m^ a = \{m h ^ + 3m^ v ) 

With the experimental values of these quantities at hand, we use a linear ansatz relating the 
three parameters (Xrhq) and the corresponding measured quantities (Y(a)). The linear approx- 
imation only holds in a limited region of the parameter space, which we estimate from earlier, 
dynamical 16 3 x 32 studies [g]. 



where the quantities Y(a) are known if we assume the lattice spacing a is known from another 
method or a-dependent if we treat a as a free parameter to be determined. Provided we are able to 
determine the J matrix and A vector, we can obtain the parameters by minimizing % 2 defined as: 



where W is the correlation matrix estimated from the measured data. We choose to use only the 
diagonal part sometimes because the data might be too noisy to give a well-behaved W. The 
quantity % 2 is a quadratic function of vector Xrhq if lattice spacing a is known and of the vector 
(moa,cp, £, a) T if a is unknown, and so it is easy to minimize analytically. J and A can be calculated 
using finite differences directly from a Cartesian set, and in order to save time we collected data 
for the minimum (seven) number of parameter sets: centered at {0.433,2.446,1.295} and with 
extent {0.1,0.1,0.02}. There is a potential problem because the RHQ parameters which we finally 
determine are actually outside of the region bounded by the 7 sets of parameters which we studied. 
However, our earlier 16 3 x 32 work suggests the region of linearity extends to this matching point. 

3. Source search and other concerns 

After studying a number of box sources we found that sources with box size 4 and 12 are the 
best to extract masses of the pseudo-scalar v\ c and vector J/y using a two state fit (t G [4, 24]), Fig. |] 
(right). However, for masses of the scalar %m and axial- vector XcU the effective mass plot, Fig. |T] 
(left) shows that without enough statistics the box size 4 source tends to give a false plateau, so we 
use the box size 12 source to determine the Xc-o an d % c \ masses via a single state fit (t € [9, 15]). 
Currently the heavy-strange data use a box source with size 8, aiming for the best plateau. The D s 




\m 1 /m 2 J 



(2.1) 



X 2 = {J-X RHQ +A-Y{a)) T W-\j-X RHQ +A-Y(a)) 



(2.2) 
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Figure 1: Sample effective mass plots of Xci (left) an d "He (right)- Red triangles: 4 sources per config (s/c) 
and box size (bs) 4; blue circles: 4 s/c, bs=12; black squares: 1 s/c, bs=4 and magenta crosses: 1 s/c, bs=12. 



andD* states are using fitting ranges t G [6,32] and [12,32] respectively. For the mass ratio m\/rri2, 
the momentum dependence is studied for both the r\ c and J/\j/ mesons and the results are quite 
consistent. We use results from the r\ c momentum dependence with the three lowest momenta. 

Other concerns such as quark propagator inversion precision are studied in detail for heavy 
quarks and the relative error for every time slice is controlled to less then 10~ 4 when a source 
placed at time zero, while light propagators are well-understood from previous studies. 



4. Analysis and results 



Let's list explicitly all the quantities used here for fitting, {X)j{m^ c + 3mj^) (2) nij/y — m^ 
(3) % (4) \{m XcQ + 3m Xcl ) (5) m Xcl - m XcQ (6) \{m Ds + 3m D *) (7) m D > - m Ds 

4.1 Heavy-heavy sector 

Using only results from the heavy-heavy states, i.e., from quantities (l)-(3), the matched RHQ 
parameters and the corresponding chiral extrapolation are shown in the table below: 





mq a 


cp 


c 


0.005 


0.410(8) 


2.356(16) 


1.270(7) 


0.01 


0.398(8) 


2.323(15) 


1.269(9) 


0.02 


0.371(9) 


2.263(14) 


1.272(9) 


0.03 


0.341(7) 


2.170(14) 


1.263(8) 


-0.00315 


0.434(9) 


2.422(18) 


1.273(9) 





niQa 


cp 


C 


0.005 


0.228(9) 


2.029(15) 


1.238(8) 


0.01 


0.217(8) 


1.998(15) 


1.237(9) 


0.02 


0.190(9) 


1.940(13) 


1.240(10) 


0.03 


0.162(8) 


1.853(13) 


1.230(8) 


-0.00315 


0.251(9) 


2.091(17) 


1.242(10) 



Table 1: The inverse lattice spacing is assumed to be 1.62 GeV (from the static quark potential with ro = 
0.50/m) for the left table and 1.73 GeV (from Q. baryon) for the right one. 

From quantities (l)-(5), we can determine the RHQ parameters as well as the lattice spacing. 



Since the states XcO and % c \ are a lot nosier the correlation matrix W (in Eq. 2.2) is not well mea- 
sured. So instead we use only the diagonal part of the correlation matrix. The chiral extrapolation 
gives aT l = 1.74(4) GeV, as shown in Fig. |2| The inconsistency between this and the result from 
the static quark potential suggests that ro is inaccurate. 
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We can make predictions for XcO an d Xci states by using J and A calculated from the mea- 
sured data and the RHQ parameters determined from quantities (l)-(3). See Fig. ||; all errors are 
propagated using the jackknife technique. When extrapolated to the chiral limit, our predictions 
are quite consistent with the experimental values, and the errors are less than one percent. This is 
encouraging and suggests that we may apply this method to make accurate predictions for other 
charmed mesons. 

3.5 | , , . . n 3.6 | , , . . n 




3.4 

-0.01 0.01 0.02 0.03 0.04 3 - 4 -ti.01 0.01 0.02 0.03 0.04 

m sea m sea 

Figure 4: The prediction of m Xc0 (left) and «», (right) in the chiral limit, with parameters determined above 
assuming aT l = 1 .73 GeV, the black lines stand for experimental values. 



Some consistency checks have been carried out for the fitting procedure, especially for the 
fitting which determines the lattice spacing. We treat lattice spacing as an input parameter, and fit 
the RHQ parameters with predictions of XcQ an d Xci extrapolated to the chiral limit as a function 
of a. Then an uncorrected, naive Xpred * s defined from: 



2 y ^Xci Ac, , ,4 y. 

Apred Z-i _n , prech v ' 

i=0 ,i 2 {m' Xa ) 

Here m p %J s means the experimental value for the %% cj meson. A plot showing the resulting %^ red 
versus the inverse lattice spacing is plotted in Fig. [| It shows good consistency that the y} pred 



p red m P h ys\i 



pn\s \ 
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minimum occurs when a is around 1.72 GeV. For all fitting procedures, such as mass fitting and 
momentum dependence fitting, we use an uncorrected fit. 

4.2 Heavy-strange sector 

The lattice spacing determined from chiral extrapolation above using heavy-heavy states is 
consistent with that determined from Q. baryon: = 1.73(2) GeV, but the errors are larger 
because of the noisy results for XcO an d Xc\ states. So we proceed to include the heavy-strange 
sector, and include the quantities (6) and (7) in the analysis to replace (4) and (5) since (6) and 
(7) are more accurately determined. The physical strange quark mass we are using is m„ = 0.036 
in lattice units, Ref. Ji2|]; but as the 24 3 data suggests a slightly different m s , this may introduce 
some systematic error. We are now studying the more accurate value m s = 0.034, so we can 
extrapolate/interpolate to the right m s assuming the dependence on the strange mass is linear. The 
results of the fitted and chiral extrapolated RHQ parameters and a with a full correlation matrix are 
listed below in Tab. ||, and the extrapolations of cT x to the chiral limit are plotted for both cases 
with full correlation matrix and diagonal correlation matrix in Fig. ||| The corresponding results 
are 1.749(14) GeV and 1.730(23) GeV respectively. 







cp 


c 


a" 1 (GeV) 


0.005 


0.241(21) 


2.052(32) 


1.240(8) 


1.722(11) 


0.01 


0.243(38) 


2.049(57) 


1.242(9) 


1.713(20) 


0.02 


0.271(30) 


2.084(42) 


1.254(9) 


1.679(14) 


0.03 


0.297(27) 


2.092(43) 


1.254(8) 


1.646(13) 


-0.00315 


0.220(28) 


2.037(42) 


1.236(9) 


1.749(14) 



Table 2: The RHQ parameters and lattice spacing determined from quantities (1)(2)(3)(6)(7), and extrapo- 
lated to the chiral limit, with % 2 from a full correlation matrix W. 

1.8 

CO 

1.75 

1.7 
1.65 

1.6 

-0.01 0.01 O02 O03 O04 -0.01 0.01 0.02 0.03 0.04 

sea sea 

Figure 5: Chiral exttapolation of inverse lattice spacing, determined from quantities (1)(2)(3)(6) and (7), 
with W a full (left) or a diagonal (right) correlation matrix. 

If we fix aT y to be 1.62 GeV then using (1)(2)(3)(6)(7) to determine the RHQ parameters will 
result in a huge % 2 jd.o.f = 146/2, which tells us the fitting fails if the wrong lattice spacing is 
used. If we set cT x to be 1.73 GeV, then x 2 jd.o.f = 1.19/2, which confirms again our observation. 
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5. Conclusion 

We have applied the RHQ action to the charmed system, both heavy-heavy and heavy-strange, 
and demonstrated that the parameters in the RHQ action can be determined with sub-percent preci- 
sion by matching several quantities to their experimental values. We discovered our lattice spacing 
from static quark potential with r<) = 0.5 fm was too large. Taking a as a free parameter we were 
then able to determine it with a few percent error. The result is quite consistent with that implied 
by the £2 baryon. In the heavy-heavy system, the XcO an d Xd states are not as well-determined 
as the J/lf/ and tj c states. We choose to use a diagonal correlation matrix when doing the four free 
parameters (three RHQ parameters and a) fitting in that case. The bare strange quark mass we are 
using in the heavy-strange run is m s = 0.036 in lattice units, which is slightly above the real one, so 
there might be a small amount of systematic error introduced into the results involving these states. 
Another heavy-strange run with different m s is underway. In conclusion, we view the application 
to the charmed system a success. And we will likely apply this method to charm-light states and 
perhaps to bottom quarks as well to explore more interesting topics. 
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